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ABSTRACT 

We make a critical assessment of models that attribute the recently detected near- 
infrared background "excess" (NIRBE) to the rcdshifted light from Population III 
objects. To supply the required 25 keV/baryon at redshift 9, Pop III massive stars 
must form with an efficiency exceeding 30% in all "minihalos" with virial tempera- 
tures above a few hundred kelvins: to avoid excessive metal pollution, most of the 
baryons once in Pop III stars must end up in intermediate-mass black holes (IMBHs). 
Gas accretion onto such IMBHs must either be inhibited or lead to early miniquasars 
with steep UV/X-ray spectra, in order not to overproduce the present-day unresolved 
soft X-ray background. In the latter case (NIRBE dominated by "X-ray quiet mini- 
quasars"), the total mass density of IMBHs at z ~ 9 must be > 50 times higher than 
the mass density of supermassive black holes observed today in the nuclei of galaxies. 
A stellar-dominated NIRBE is less economical energetically: > 5% of all baryons in 
the universe must be processed into Population III stars. We survey various aspects 
of the Population III hypothesis for the NIRBE, and show that the ionizing photon 
budget required to account for the NIRBE is much larger than that required to explain 
the high electron scattering optical depth measured by the WMAP satellite. 

Key words: cosmology: diffuse radiation - galaxies: formation - intergalactic medium 



1 INTRODUCTION 

Independent measurements based on COBE/D1KBE and 
JiiTS/NIRS data provide evidence for an excess isotropic 
emission of extragalactic origin in the wavelength band 1.2 
to 4 ^m (Dwek & Arendt 1998; Gorjian et al. 2000; Wright 
& Reese 2000; Wright 2001; Cambresy et al. 2001; Mat- 
sumoto et al. 2004). The inferred flux appears too bright to 
be accounted for by the integrated light from faint galaxies, 
and shows a spectral discontinuity from the resolved opti- 
cal background at ~ 1/im (Wright & Reese 2000; Cambresy 
et al. 2001; Matsumoto et al. 2004). A similar conclusion 
is obtained from the analysis of near-infrared background 
(NIRB) fluctuations on angular scales ranging from half a 
degree (DIRBE) to sub-arcmin (2MASS), namely that nor- 
mal galaxies cannot account for the observed angular power 
spectrum (Kashlinsky & Odenwald 2000; Kashlinsky et al. 
2002). It has been sugggested that the observed NIRB spec- 
tral excess and its brightness fluctuations may be the signa- 
ture of the redshifted UV light from Population III (zero- 
metallicity) star formation at z ~ 10 (Santos et al. 2002; 
Salvaterra & Ferrara 2003; Magliocchetti et al. 2003; Kash- 
linsky et al. 2004; Cooray et al. 2004). Alternative possibil- 
ities involve Eddington-limited emission from an early pop- 



ulation of accreting miniquasars (Cooray & Yoshida 2004), 
or, perhaps less exotically, flux missed by standard photom- 
etry in the outer, lower surface brightness parts of galaxies 
(Bernstein et al. 2002). 

In this paper we assess the merits of Pop III objects 
as the root cause of the NIRB excess (NIRBE). We show 
that the energetic requirements are uncomfortably high, and 
discuss the astrophysical implications of an early epoch of 
massive Pop III star formation. 



2 A POP III MODEL FOR THE NIRBE 
2.1 Simple energetics 

In the near IR, out of ~ 35 nWm~ 2 sr _1 of integrated flux 
from 1 to 4 /im, about Ij — 13.5±4.2 nWm~ 2 sr _1 are con- 
tributed by the Cambresy et al. (2001) DIRBE measurement 
in the J-band at 1.25 /xm. 1 This number is quite uncertain, as 
it depends on zodiacal light subtraction (Kelsall et al. 1998; 

1 Throughout this paper we will quote surface brightness values 
integrated over the bandwidth, rather than the more common 
brightness per logarithmic bandwidth. 
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Wright 1998). Wright's zodiacal light model yields a lower 
value for the J-band NIRB, Ij = 7.2 ± 4.0 nWrn^sr" 1 
(Wright 2001). The contribution from known sources ob- 
tained by integrating the deepest galaxy counts at 1.1 [im 
amounts to 2.6t°\l nWm~ 2 sr" 1 (Madau & Pozzetti 2000). 
The leftover light that is often attributed to Pop III stars 
then ranges from a few to about 10 nWm _2 sr _1 . To be 
conservative, we shall adopt in the following a representa- 
tive value for the NIRBE in the J-band of 2.5 nWm -2 sr _1 . 
This corresponds to an observed excess energy density of 

u j = 10~ 15 ergs cm -3 . (1) 

Assume, for simplicity, that all this energy was emitted as 
Lya radiation 2 by very massive Pop III stars with a narrow 
redshift distribution centered around (z) = 9. This is the 
redshift at which the data requires a rapid transition from 
Pop III to normal Population II stars (Salvaterra & Ferrara 
2003): a spectral jump in the NIRB then occurs at the red- 
shifted wavelength (l + (z)) 1216 A. The energetics would be 
even more severe if the sources responsible for the NIRBE 
were broadly distributed in the range z — 9 — 20 (say), or 
if the initial stellar mass function (IMF) was not dominated 
by very massive stars. The comoving Lya energy density 
radiated at (z) = 9 is then 



U a = (l + (z))uj = 10 uj, 



(2) 



corresponding to 25 keV per cosmic baryon. Zero-metallicity 
stars in the range 80 < m* < 500 M emit AT 7 = 70, 000 - 
90, 000 photons above 1 ryd per stellar baryon over a lifetime 
of t ms = 2 - 3 x 10 6 yr (Schaerer 2002). In the following, 
we will adopt a mid-range value of iV 7 = 80, 000 ionizing 
photons per stellar baryon. In the Case B approximation, a 
fraction 0.7 of these is converted into Lya . The total mass 
density of baryons that needs to be processed by (z) = 9 
into Pop III stars to produce the NIRBE is then simply 
estimated as 



p* = 



m p U a 
0.7N-,E a 



2.7 x 10 s M Mpc" 



(3) 



where m p is the proton mass and E a — 10.2 eV is the en- 
ergy corresponding to the Lya transition. Note that this 
mass density depends only weakly on the IMF of massive 
(> 80 Mq) Pop III stars. Dividing p, by the critical density 
one gets 3 



n* = 0.002 = 0.045 Q b 



(4) 



i.e. at least 5% of all baryons in the universe must be con- 
verted into Pop III stars (see also Santos et al. 2002). This 
is energetically and astrophysically daunting. It means Pop 
III star formation at early times must be comparable with 
the mass density in stars observed today, (2.3 ± 0.3) x 10~ 3 
(for a Kennicutt IMF; Cole et al. 2001). 

2 This is energetically advantageous, since most of the power of 
massive Pop III stars is emitted above 1 ryd, and in Case B re- 
combination approximately 50% of this energy is re-emitted as 
Lya. 

3 Throughout this paper we adopt a flat ACDM back- 
ground cosmology with parameters (Oa, fijf, Hj,, n, as, h) = 
(0.3,0.7,0.045, 1,0.9,0.7). 



2.2 Star formation efficiency 

The minimum mass scale for the gravitational aggregation 
of cold dark matter particles is negligibly small. In the Press- 
Schechter formalism, the mass fraction in dark matter halos 
above mass M at redshift z can be written as 



F(> M\z) = erfc 



S c (z) 



\/2a M {z) 



(5) 



where om(z) is the linear theory rms density fluctuation 
smoothed with a 'top-hat' filter of mass M at redshift z, and 
S c (z) is the critical threshold on the linear overdensity for 
spherical collapse at that redshift (Press & Schechter 1974). 
Baryons, however, respond to pressure gradients and need 
to contract and cool in order to fragment into stars. It is 
useful here to identify two mass scales: (1) a molecular cool- 
ing mass Mh 2 above which gas can cool via roto-vibrational 
levels of H 2 and contract, Mh 2 « 10 s [(1 + z)/10p 3/2 Mq 
(virial temperature in excess of 200 K; Machacek et al. 2003); 
and (2) an atomic cooling mass Mh above which gas can 
cool efficiently and fragment via excitation of hydrogen Lya, 
M H « 10 s [(1 + z)/10]~ 3/2 Mq (virial temperature above 
10 4 K). Figure shows the fraction of the total mass in the 
universe that is in collapsed dark matter halos with masses 
greater than Mh 2 and Mh at the epochs of interest here. 
Let us define a star formation efficiency /* as the fraction of 
collapsed baryons that cool and form Pop III stars, 



n b gF(> Mh 2 |9)' 



(6) 



Here g < 1 is a correction factor accounting for the fact that 
the gas fraction (gas mass/virial mass) in "minihalos" con- 
densing due to H2 cooling may be lower than the universal 
cosmic value Qb/flm- In the absence of a UV photodissoci- 
ating flux and of ionizing X-ray radiation, three-dimensional 
simulations of early structure formation yield gas fractions 
g = 0.6 — 0.7 (Machacek et al. 2003). The same simulations 
also show that the amount of cold, dense gas available for 
star formation exceeds 20% only for halos above 10 6 Mq. 

From equation 01 and Figure Q it is obvious that if the 
whole J-band excess is caused by Pop III stars, then one 
needs to invoke high star formation efficiencies, /, w 0.3, 
in all halos above the molecular cooling mass. Pop III star 
formation that was limited only to halos above the atomic 
cooling mass could make a significant contribution to the 
NIRBE only for /, ~ 1, which is implausible. Our estimate 
of /, agrees with that of Santos et al. (2002) but is higher 
than the value of 10% derived by Salvaterra & Ferrara (2003) 
for a very top-heavy IMF. While comparably high efficien- 
cies appear to be needed for globular clusters to survive 
gas expulsion and remain gravitationally bound (Goodwin 
1997), it is not clear how such large values may be reached 
in metal-free gas restricted to relatively inefficient H2 cool- 
ing. It is also surprising that negative radiative, mechanical, 
and chemical "feedback" effects would not be choking off 
the formation of such a large number of Pop III stars, given 
the fragile nature of these molecules (Haiman et al. 2000; 
Ricotti et al. 2002; Scannapieco et al. 2002). 
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Figure 1. Solid lines: Total mass fraction in all collapsed dark 
matter halos above the molecular cooling and the atomic cool- 
ing masses, Mjj 2 and Mh, as a function of redshift. Dashed line: 
NIRB (nWm _2 sr _1 ) due to redshifted Lya emission from Popu- 
lation III stars at redshift z (for /* = 0.3, g = 0.65, = 80, 000, 
and F(> Mh 2 \z), see text for details). 

2.3 Reprocessing of ionizing radiation 

The interpretation of the J-band NIRBE as redshifted 
Lya from (z) = 9 requires the emission of U a /(0-7E a rib) ~ 
3, 500 Lyman continum photons per baryon in the uni- 
verse. Here rib = 2.5 x 10~ 7 cm~ 3 is the mean cosmologi- 
cal baryon density today. The assumption that most of this 
radiation gets absorbed and converted into Lya within the 
host halo then implies that every baryon in the halo that 
was not converted into stars must recombine on average 
N ICC = iV 7 /*/(l - /») « 35,000 times over the lifetime 
of a Pop III star cluster. The volume- averaged hydrogen re- 
combination timescale for halo ionized gas at temperature 
10 4 K with mean overdensity 8 = 200 feoo relative to the 
background baryon density is 

t ICC « 3 x 10 6 yr (i±£) ' S^ (7) 

where Ch = { n p)/{ n p} 2 takes into account the degree of 
dumpiness of photoionized halo gas (here (n p ) is the mean 
number density of ionized hydrogen in the halo). Then 
from iV rcc = (tms/trcc) « [(1 + z)/10] 3 fcoo Ch , one derives 
a clumping factor Ch ~ N tec , i.e. the Lyman continuum 
photons must be absorbed in the densest molecular clouds 
around the sites of star formation (see also Santos et al. 
2002). A small ionizing-photon escape fraction, / osc , from 
minihalos into the IGM is in contrast with the recent re- 
sults of 1-d radiation hydrodynamical calculations of the 
evolution of H n regions around Pop III stars (Whalen et al. 
2004). These show that the ionization fronts exit the halo 
on timescales much shorter than the stars' main sequence 
lifetimes: the photoevaporative flows leave dark halos with 
a low gas content and produce escape fractions / csc > 0.95. 
If this were the case, reprocessing of ionizing radiation into 
Lya would still occur in the low-density IGM (Santos et al. 



2002; Haardt & Madau 1996). Note that Case B recombina- 
tion predicts approximately 10% of the Lya luminosity to 
be reprocessed into redshifted Ha line emission at 6.5 /im. 



2.4 Metal enrichment 

Non-rotating Pop III stars in the mass window 140 < 
m* < 260 Mq will disappear as pair-instability supernovae 
(PISNe; Bond et al. 1984), leaving no compact remnants 
and polluting the universe with heavy elements. Stars with 
40 < m* < 140 M and m* > 260 M Q are predicted instead 
to collapse to IMBHs: the mass of metals (mostly oxygen) 
ejected by PISNe is approximately equal to m*/2 (Heger 
& Woosley 2002). Therefore if stars giving birth to PISNe 
were responsible for the NIRBE, they would overenrich the 
universe to a mean metallicity, (Z) — Qz/^t = 0.5fi*/fib « 
0.025, in excess of solar already by redshift 9 or so. Even 
allowing for only partial mixing of the ejecta, it is difficult 
to see how all these metals could remain undetected: the 
unusual nucleosynthetic signature of PISNe would appear 
in the metal abundances of second-generation stars, 4 while 
absorption studies of the Lya forest at z ~ 3 would reveal 
cosmic metallicities far in excess (by > 2 orders of magni- 
tude) of the observed values (Schaye et al. 2003). A similar 
argument has recently been used to strongly constrain the 
possible contribution of PISNe to the reionisation of the uni- 
verse (Daigne et al. 2004). 



2.5 Miniquasars 

As already pointed out by Santos et al. (2002), if distant 
Pop III stars make a significant contribution to the NIRBE, 
then perhaps the most conservative hypothesis is to assume 
that most of the baryons once in Pop III stars ended up in 
IMBHs. If these seed holes were able to increase their initial 
mass via gas accretion and shine as miniquasars, and a sig- 
nificant fraction of the radiated energy was emitted above 
and close to the hydrogen Lyman edge, then they - and not 
their metal-free stellar progenitors - would dominate the UV 
background flux (Madau et al. 2004). Let us denote with /uv 
the fraction of the bolometric power radiated by miniquasars 
that is emitted as hydrogen-ionizing radiation with mean 
energy (E). Then the number of Lyman continuum pho- 
tons produced per accreted baryon is iV 7 = rjm p c 2 /uv / (E) , 
where r\ is the radiative efficiency. Accretion via a thin disk 
predicts most IMBHs at early times to be rapidly rotat- 
ing with radiative efficiencies exceeding 10-15% (Volonteri et 
al. 2005). If the shape of the emitted spectrum from mini- 
quasars followed the mean spectral energy distribution of 



4 It has been suggested that the relative abundances in extremely 
metal-poor Galactic halo stars may show evidence for PISN en- 
richment below metallicities [Fc/H]= —3 (Qian & Wasserburg 
2002). The inferred fraction of baryons in the universe processed 
into PISN progenitors is, however, only < 10 — 4 (Oh et al. 2001), 
approximately 500 times lower than the value needed to explain 
the NIRBE (see eq.gj. 
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the quasar sample in Elvis et al. (1994), then /uv = 0.3 
and (E) = 3ryd. One then estimates iV 7 « 10 6 T7o.i5, more 
than 10 times larger than the typical value for Pop III stars. 
Under the assumption that all these photons were converted 
into Lycv radiation at (z) = 9, the total mass density that 
must be accreted onto IMBHs to explain the NIRBE is 



Pace 



m p U a 
0.7 'N-yE a 



2 x 10 7 i]o l 5 M Mpc" 



(8) 



In the case of accretion at the Eddington rate, the hole 
mass exponentiates over a Salpeter timescale, ts = (4.5 x 
10 s yr) 77/(1 — 77). By redshift 9, the mass density of Pop 
III holes could increase by as many as 6 e-foldings (taking 
r\ = 0.15), or about a factor 400. A miniquasar model for 
the NIRBE may then require an initial IMBH progenitor 
mass density of only p* « p acc /400 w 50, 000 M Mpc~ 3 , 
> 5000 times lower than that given in equation @. Note 
that our results differ from those presented by Cooray & 
Yoshida (2004). These authors fix the amount of Pop III 
star formation occurring at early times and then use the 
observed NIRBE to constrain the fraction of stellar mass 
converted into Eddington-limited seed IMBHs to be < 10%. 
To avoid overenrichment we take the seed hole fraction to be 
unity instead, and allow the mass density of Pop III holes to 
exponentiate via gas accretion, thus decreasing the amount 
of baryons that needs to be converted initially into Pop III 
stars. 

The scenario described above is, in principle, an attrac- 
tive solution to the NIRBE problem, as it is more economical 
energetically and does not demand very high star formation 
efficiencies down to minihalo masses. Still, plenty of cold 
gas must be made available at early times for IMBHs to ac- 
crete, driving the total mass density of IMBHs at 2 ~ 9 to 
more than 50 times higher than the inferred mass density of 
supermassive black holes today (Yu & Tremaine 2002). The 
main issue, however, is that miniquasars powered by IMBHs 
are expected to be hard X-ray emitters. Hard X-ray photons 
produced in the (0.5 — 2.0)(1 + (z))keV band by our pop- 
ulation of distant miniquasars would redshift without ab- 
sorption and saturate the unresolved soft X-ray background 
(SXB). For a simple power-law Fe oc E~ a with a > 0, the 
relative contribution to the NIRBE and the SXB is 



TtSXB 



0.5 



(4M( 



5keV 
,13.6eV 



/ 20keV \ ] 
U3.6eV/ 



.(9) 



The energy density associated with the unaccounted com- 
ponent of the present-day 0.5 — 2keV SXB does not ex- 
ceed «sxb ~ 2 x 10 -18 ergscm -3 (Dijkstra et al. 2004), 
500 times smaller than the NIRBE in the J-band. Equa- 
tion is consistent with the observations only for power- 
laws steeper than a > 2.3 (see also Salvaterra et al. 2005). 
While the photon energy distribution of putative high-z 
miniquasars is very uncertain, the spectra of 'ultraluminous' 
X-ray sources in nearby galaxies appear to require both a 
soft component (well fit by a cool multicolor disk black- 
body with fcT max ~ 0.15 keV, which may indicate IMBHs; 



5 This has the same 2500 A— 2 keV spectral slope, tiox 
the Sazonov et al. (2004) template. 



1.4, of 



Miller & Colbert 2004) and a non-thermal power-law com- 
ponent of comparable luminosity and slope a ~ 1. Even if 
miniquasars powered by accreting IMBHs were so "X-ray 
quiet" as to make a dominant contribution to the NIRBE 
while satisfying the SXB constraint, only a small fraction 
f w of the energy released could be used to drive an out- 
flow and be ultimately deposited as thermal energy into the 
IGM. Assuming rapid, homogeneous "preheating" at red- 
shift 9 (Benson & Madau 2003), the comoving thermal en- 
ergy density introduced in the IGM would be Uigm = fwU a - 
At these early epochs, inverse Compton scattering would 
transfer all the energy released to the cosmic microwave 
background (CMB), producing a y-distortion to its spec- 
trum, y = [/igm/(4[/ C mb) = 6 x 10~ 4 f w . The COBE satel- 
lite measured y < 1.5 x 10~ 5 (Fixsen et al. 1996), implying 
f w < 0.025. 



2.6 IMBHs in galaxy halos 

We are therefore led to consider the possibility of a large 
population of IMBHs "wandering" in galaxy halos and un- 
able to accrete significant amounts of material over a Hubble 
time. As the NIRBE would now be dominated by their stel- 
lar progenitors, the total mass density of remnant IMBHs, 
« 0.05f2f,, would exceed the mass density of the supermas- 
sive variety found today in the nuclei of most nearby galax- 
ies, fisMBH = (2.1 ±0.3) x 10~ 6 (Yu & Tremaine 2002), by 3 
orders of magnitude. The abundance of IMBHs - assuming 
they behave similarly to collisionless dark matter particles 
- is a function of the environment, being higher (positively 
biased) within a massive galaxy bulge that collapsed at early 
times, and lower (antibiased) within a dwarf galaxy collaps- 
ing at the present-epoch. This effect can be quite easily 
quantified within the extended Press-Schechter formalism, 
and is expected to be weak for progenitor halos of mass Mh2 
collapsing at 2 = 9 and merging at later times into a "Milky 
Way" halo (Madau & Rees 2001). As the host of IMBHs ag- 
gregate hierarchically into more massive systems, dynamical 
friction against the visible and dark matter background will 
favor accretion of larger mass subunits and infalling satel- 
lites with low initial angular momenta: their cores will merge 
and undergo violent relaxation. Most IMBHs will not sink to 
the center and will be left wandering in galaxy halos (Madau 
& Rees 2001; Volonteri et al. 2003; Islam et al. 2003), where 
they can most easily escape detection. 



3 DISCUSSION 

The Pop III hypothesis for the NIRBE has appeared of- 
ten in studies of the extragalactic background light, the first 
stars and their importance for reionization, the duration and 
extent of metal- free star formation (Kashlinsky 2004). It is 
only fair to point out that the very existence of this excess is 
open to question, because of the difficulty of accurately sub- 
tracting the contribution of bright foregrounds. After sur- 
veying many aspects of the Pop III model for the NIRBE, it 
is hard to reach any firm conclusion regarding its tenability, 
the unknown initial mass function of Pop III stars being one 
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of the many adjustable theoretical parameters. Most previ- 
ous studies have emphasized the possibility that both the 
NIRB fluctuations and mean levels may provide direct in- 
formation on the epoch of the first stars. While a pregalactic 
contribution to the NIRB at a level of a few nWm _2 sr~ 1 
cannot be ruled out, we offer here an assessment of the main 
weaknesses of the Pop III hypothesis. For a given amount of 
material converted into stars or accreted onto black holes, a 
factor (1 + z) is lost to cosmic expansion in the total back- 
ground light observed at Earth when converting from ob- 
served to radiated (comoving) luminosity density. An early 
epoch of production of the NIRBE thus makes the energetic 
requirements even more extreme. A very top-heavy stellar 
IMF helps with the energy budget but overenriches the uni- 
verse with heavy elements unless most baryons once in Pop 
III stars end up in IMBHs. A fine tuning of the IMF is then 
required to avoid the PISNe mass window. Radiately effi- 
cient accretion onto such a large population of IMBHs will, 
however, easily overproduce the SXB. 

The recent detection by the Wilkinson Microwave 
Anisotropy Probe ( WMAP) satellite of a large optical depth 
to Thomson scattering, r e = 0.17±0.04, suggests that ioniz- 
ing sources were already very abundant at redshifts z > 10 
(Spergel et al. 2003). If the universe is suddenly reionized 
at redshift z r , the Thomson scattering optical depth to z T 
is T e = \l) . <jt n e (t) cdt, where or is the Thomson cross 
section and n e (t) is the proper free-electron density. The 
mean number of recombinations per hydrogen atom in the 
IGM since redshift z r is A/" rC c = J t *,° dt/t Tec . As long as 

the temperature of the IGM remains constant, T = 10 4 K, 
with redshift, J\f rac can be written as a function of r e alone 
(Dijkstra et al. 2004) as 



. _ qb(T) Cigm 

^Vrec — Te 



<7T C 



Cigm (r e /0.i 



(10) 



where Qs is the recombination coefficient. Here Cigm is the 
clumping factor of the diffuse IGM at high redshift and is 
expected to be of order unity (we purposely exclude halo gas 
in our definition of Cigm , as the consumption of ionizing ra- 
diation by halo gas is accounted for by the escape fraction 
/esc). To keep intergalactic gas ionized, one then needs the 
production of A/" 7 = (1 + A/" re c)//esc Lyman continuum pho- 
tons per atom in the IGM. A J-band excess at a level of a few 
nWm -2 sr -1 , interpreted as redshifted Lycv from redshift 
10, implies the emission of approximately 4,500 H-ionizing 
photons per atom in the universe. Even for ionizing-photon 
escape fractions of order a percent or so, the ionizing photon 
budget required to account for the NIRBE is clearly much 
larger than that required to explain the WMAP results. 
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